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Abstract: We report radio continuum observations with the Australia Telescope Compact 
Array of two molecular clouds. The impetus for these observations is a search for syn- 
chrotron radiation by cosmic ray secondary electrons/positrons in a region of enhanced 
density and possibly high magnetic field. We present modelling which shows that there 
should be an appreciable flux of synchrotron above the more difl*use, galactic synchrotron 
background. 

The starless core G333. 125-0.562 and infrared source IRAS 15596-5301 were observed at 
1384 and 2368 MHz. For G333. 125-0.562, we flnd no signiflcant levels of radio emission 
from this source at either frequency, nor any appreciable polarisation: we place an upper 
limit on the radio continuum flux from this source of 0.5 mJy beam~^ at both 1384 and 
2368 MHz. Due to the higher than expected flux density limits, we also obtained archival 
ATCA data at 8640 MHz for this cloud and place an upper limit on the flux density of 
50 /iJy beam~^. Assuming the cosmic ray spectrum is similar to that near the Sun, and 
given the cloud's molecular density and mass, we place an upper limit on the magnetic 
field of 500 /xG. 

IRAS 15596-5301, with an RMS of 50 /iJy beam~^ at 1384 MHz, shows an HII region 
consistent with optically thin free-free emission already detected at 4800 MHz. We use 
the same prescription as G333 to constrain the magnetic field from this cloud to be less 
than 500 /iG. We find that these values are not inconsistent with the view that magnetic 
field values scale with the average density of the molecular cloud. 

Keywords: synchrotron radiation, cosmic rays, molecular clouds, starless cores: general — molecular 
clouds: individual: G333.125-0.562, IRAS15596-5301 



1 Introduction 

The Galactic synchrotron emission is usually thought 
to be due entirely to primary cosmic ray electrons and 
electron/interstellar medium (ISM) interactions. How- 
ever, a significant contribution could also be made by 
electrons and positrons (hereafter e^) produced as sec- 
ondaries by interactions of cosmic ray (CR) nuclei with 
atomic nuclei, in molecular clouds where conditions 
(such as magnetic field strength and density) could be 
favourable for the production of synchrotron photons 
from the secondary e^. 

CRs play an important role in molecular cloud evo- 
lution by partially ionizing even the cold molecular gas, 
and thereby affecting the dynamics of cloud collapse 
by coupling the magnetic field to the partially ionized 
gas. This in turn could result in amplification of the 
ambient magnetic field during cloud collapse, and give 
rise to a correlation between d ensity of molecular gas 
and magnetic field. Crutcher ( p,999J has found just 
such a relation between magnetic field and average H2 
number density, (cm~^), of individual molecular 
clouds. 



We recently conducted a radio continuum study 
with the Australia Telescope Compact Array (ATCA rj 
to search for evidence of synchrotron emission due to 
secondary in molecular clouds. As will be shown, 
the expected emission is linearly proportional to tar- 
get mass, and should have a negative spectral index 
a ^ —0.7 (for Si^ oc z/", where a is the spectral in- 
dex) . We therefore selected massive clouds such as the 
Galactic Centre (GC) giant molecular cloud (CMC), 
Sagittarius B2 (Sgr B2) because of the higher expected 
CR density in the region. Because of the consider- 
ably more complex structure of this cloud, the results 
of the observations of Sgr B2 will be presented else- 
where (Jones et al., in preparation; see also Crocker, 
et al. 2007). We also selected nearby 3-4 kpc) 
cores of giant molecular cloud (CMC) complexes, and 
in order to maximise the probability of detection, we 
selected clouds which have little or no known ther- 
mal radio emission. The sources selected were the in- 
frared source, IRAS 15596-5301 (hereafter IRAS 15596) 
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■"^The Australia Telescope Compact Array is a part of the 
Australia Telescope which is funded by the Commonwealth 
of Australia for operation as a National Facility managed 
by CSIRO 



2 



Publications of the Astronomical Society of Australia 



( |Garay et al.|[2004| , for its low levels of radio con 
tinuum emission at 4800 MHz and no spectral infor 
mation, an d G333. 125-0.562 (hereafter G333) ( |Garay 
et al.|20()2| , because of its status as a starless core and 
lack of radio continuum observations. 



Table 1: Distances to and physical parameters for 
the two clouds. 



IRAS15596' 



Mass (M0) 
uh^ (cm-^) 
T(K) 
d (kpc) 
R{pc) 



2.3x10^'^ 
2 X 10^ 
13^ 
3.6 
0.34 



1.4x10^ 
4 X 10^ 
- 27 
4.6 
0.47 



Garay et al. (2004); 



et al. 



(2007) 



Garay et al. (2002); ^ Lo 



The densities of the two clouds, are listed in Ta- 
ble [l] and using the Crutcher prescription, we expect 
magnetic field values of 400t^^° and 550t?°° fiG for 
G333 and IRAS 15596 respectively. The error for the 
expected values is derived from the diQ.QS u ncertainty 
of the exponent quoted by I Crutcher (1999). 



The production rate of secondary electrons and 
positrons depends only on the spectrum and intensity 
of cosmic ray nuclei, and the density of the interstellar 
matter. We use the production rate of electrons and 
positrons qm{E), per solar mass of interstellar matter 
per unit energy (Mq^ GeV~^ for the cosmic ray 

spectrum and composition observed above the Earth's 
atmosphere based on Figure 4 of |Moskalenko Strong] 
(|1998| (we note here that the choice of CR spectral 
index and normalisation is very important, since the 
clouds are so close, there is no reason to suppose that 
the spectrum should deviate significantly from that ob- 
served above our atmosphere). The production spec- 
trum of electrons and positrons per unit volume per 
unit energy at position r is then the product of the 
density of interstellar gas at position r multiplied by 
the production rate per unit mass (e^ cm~^ GeV~^ 

q{E,r) =nH^{Y)qm{E)2mH. (2) 
For moderate molecular cloud densities uh > 10^ cm" 
and magnetic fields B > 10~^ G, the relatively short 
energy loss time justifies neglecting diffusive transport 
of electrons within the cloud. Additionally, we note 
that the small radii of the clouds themselves, should 
not prevent the ^ GeV CR protons from freely enter- 
ing the clouds. Hence we assume that the clouds are 
'bathed' in GeV CRs, and the charged secondary par- 
ticles produced decay and cool within the confines of 
the cloud. 

Having obtained the production spectrum, one then 
readily obtains by numerical integration the ambient 
number density of electrons and positrons, per unit en- 
ergy, n^(^,r)(e^ cm~^ GeV~^), at various positions 
r within the molecular cloud complex: 



where dE/dt is the total rate of energy loss of electrons 
at energy E due to ionization, bremmstrahlung and 
synchrotron emission (because of the energies involved, 
we neglect positron annihilation, and assume electrons 
and positrons suffer identical energy losses) . Electrons 
lose energy by ionization losses in neutral molecular 
hydrogen at a rate (in GeV s~"^) 



n{E,r) 



J^q{E,r)dE 
dE/dt 



(3) 



dE 

dt loniz 

dE 

dt bremss 

dE 

dt synch 



-5.5 X 10-1^ f-^^'jxani+e.SS), (4) 
\ 1 cm / 

-■■=-"'-'(r<lv)-(il^) 



= -1.0 X 10" 



j^) X7^ (6) 
1 gauss / 

where B± is the component of magnetic field perpen- 
dicular to the electron's direction. For an isotropic 
electron population {B±) = 0.78B. The synchrotron 
emissivity, ji^ (in erg cm~^ s""*^ sr""*^ Hz""*^), is then 
calculated using stand ard formulae in synch rotron ra- 
diation theory ( ^Rybicki Lightman TT979t 



AnnieC^ 



(7) 



4.19 X 10^(^/mec^)^(- 



1 Gauss 



: 4.8 X 10" 



niec 8.18 X 10"'^ 



Hz, (8) 

esu, (9) 
erg,, (10) 



F{x) 



K^(i)di. (11) 



where ix) is the modified Bessel function of order 

3 

5/3. 




8 9 10 

log[v/Hz] 



Figure 1: Log- log plot of the expected flux from 
the two cores chosen for the observations; solid line 
is G333, dotted line is IRAS15596. The parame- 
ters used are listed in Table 1, along with the value 
of the magnetic field as predicted by the [Crutcher 



( |1999 l) scaling relation. Additionally, the local CR 
spectrum is assumed (for reasons, see text). 
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Figure [T] shows the expected synchrotron emission 
from the two clouds using the above outUned method, 
the physical parameters listed in Table 1 , and the Crutcher 
scaling relation. This shows that a sensitivity on the 
order of 100 /xJy (at centimeter wavelengths) should 
provide a robust detection of these clouds at centime- 
ter wavelengths. 



2 Observations 

2.1 Observations Data Reduction 

Single pointing ATCA radio continuum observations 
were made on 2005 December 4 — 5*^, using the 1.5C 
configuration utilising five of a possible six antennae 
and covering east- west baselines from 76.5 meters to 
1.4 km. Observations were made simultaneously at 
the frequencies 1384 MHz (A = 20 cm) and 2368 MHz 
(A = 13 cm), each spanning 128 MHz in bandwidth 
and recording all four polarisation products. At these 
frequencies, the FWHM of the ATCA beam is 34.5' 
and 20.1' at 1384 and 2368 MHz respectively. The to- 
tal integration time for each source was 720 minutes 
obtained from 50 minute scans over 2 hours on the 4*^ 
for G333. 125-0.562 and 2 hours on the 5*^ for IRAS 
15596-5302. The remaining observations were gained 
by using 15 minute cycles over 10 hours on the 5*^, in 
order to provide excellent ('iz^'u)- coverage. The calibra- 
tor PKS 1613-586 was observed before and after every 
on-source scan in order to correct the amplitude and 
phase of the interferometer data for atmospheric and 
instrumental effects. The flux density was calibrated 
by observing PKS 1934-638 (3C84), for which the flux 
densities of 14.94 and 11.59 Jy at 1384 and 2368 MHz 
respectively were adopted. Due to PKS 1934-638 be- 
ing below the horizon, the array was system calibrated 
using PKS 0823-500, and PKS 1934-638 was observed 
at the end of both of the observing runs, and the flux 
densities for both PKS 0823-500 and PKS 1613-586 
were bootstrapped from PKS 1934-638. 

Data reduction was performed using the MIRIAD 
software package using standard calibration procedures, 
except as noted above. Imaging was then performed by 
Fourier transforming the interferometer data using the 
MIRIAD task invert The images were then CLEANed 
using the task maxen, where the flux was constrained 
using a total flux estimate gained from a CLEANed 
image gained using the task clean. The synthesised 
(FWHM) beams were at 1384 and 2368 MHz respec- 
tively: 20''. 3 X 18''. 7 and 11". 3 x 9". 9 for G333 and 
19" X 10". 5 and 15" x9".5 for IRAS15596. The noise in 
the images were determined by adding in quadrature, 
the standard deviation of the flux of the non-primary- 
beam-corrected image, integrated within a number of 
regions roughly corresponding to the beam size in ar- 
eas of the image which lack an obvious source. The 
RMS noise was determined to be 0.5 mJy beam~^ at 
both 1384 and 2368 MHz for the G333 cold core, and 
30 and 50 /xJy beam~^ at 1384 and 2368 MHz respec- 
tively for IRAS15596. The RMS for the IRAS source 
are close to the theoretical noise limit of the telescope, 
whereas for the G333 cold core, it is ^ 10 times the 



theoretical noise limit due to strong emission nearby. 
2.1.1 Additional Data 

In order to better characterise any emission from the 
chosen clouds, we have obtained additional, archival 
data: 843 MHz Molongolo Observatory Synthesis Tele- 
scope (MOST) data for both G333 and IRAS15596; 
and archival 8640 MHz ATCA data for G333. 

843 MHz MOST data was obtained from the Syd- 
ney Unive rsity Molongolo Syn thesis Survey (SUMSS) 
catalogue ( |Mauch et aL]|2003| ). SUMSS has a synthe- 
sised beam of 43 x 43 sin~ (S), where S is the decli- 
nation of the source, and a limiting peak brightness of 
between 6 a nd 10 mJy beam~^ , depending on source 
declination ( |Mauch et"aLl|2003 ). 

We also obtained archival ATCA data for the G333 
cloud for two reasons: firstly to check for low level, 
optically thick emission at higher frequencies, and sec- 
ondly (as will be show in §4) the RMS noise limit for 
the lower frequencies are much higher than the ex- 
pected emission shown in Figure The data was from 
Project Code C483, taken in 199^ It used the 6A and 
C configurations of the ATCA, utilising all 6 possible 
antennae, with baselines from 153 and 337 m to 6 km, 
for the 6C and 6A configurations respectively. The ob- 
servations were centered at a frequency of 8640 MHz 
(A = 3 cm) spanning 128 MHz in bandwidth (4800 
MHz (A = 6 cm) observations were taken simultane- 
ously, however, these were in spectral line mode, so 
were not used). It is unclear how many polarisation 
products were recorded, however only Stokes I is used 
here. 

These data were reduced using the steps outlined 
above in miriad, using PKS 1934-638 as the primary 
calibrator (assuming the same flux densities noted above), 
and PKS 1600-40 as the secondary calibrator (the sec- 
ondary calibrator was chosen by plotting the {u, v)- 
distance as a function of time for the calibrators and 
source field, and selecting that which had observations 
of the caibrator interspersed). At 8640 MHz, the pri- 
mary beam is 5.5' and the resulting synthesised beam 
is 1.2" X 0.8" at a position angle of 24°. The 1-a RMS 
noise for the image was found to be ^ 50 /xJy beam~^, 
just above the expected noise level. 



3 Results 

3.1 G333.125-0.562 

Figure |2] shows the total intensity image of the region 
around the G333 cold core from the 1384 MHz ATCA 
observatio ns, overlaid with the 1.2 mm contours as de- 
scribed in Garay et al. ( |2004 ). The bright emission 
in the image is a large ridge of emission known as the 
RCW106 giant molecular cloud (CMC). It is obvious 
that even with a logarithmic transfer function applied 
to the image, there is no source coincident with the 

1.2 mm dust emission from G333 as observed by |Garay| 
|etal.||2004| down to the 1-a level at 843, 1384 or 2368 



^It can be accessed from the ATCA on-line archive: 
http://atoa.atnf.csiro.au/ 
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MHz. There was no significant polarised emission de- 
tected at either 1384 or 2368 MHz. 




Right Ascension 



Figure 2: G333.125-0.562: ATCA 1384 MHz to- 
tal intensity image overlaid with (black) 1.2 mm 
(Garay, 2004) contours at 0.18, 0.36, 0.54, 0.90, 
1.36 and 1.80 Jy beam~^. The intensity image uses 
a logarithmic transfer function - \og{v al — rain) - 
to bring out emission on all scales. The beam for 
the 1384 MHz emission is located in the lower left- 
hand corner. The G333 cold core is labelled. 

Since the noise level from the 1384 and 2368 MHz 
images does not match the expected source strength, 
due to the extended emission in the beam, we imaged 
the region using the archival 8640 MHz ATCA data 
described above. The 8640 MHz data, whilst at a fre- 
quency where it could be expected to be confused with 
thermal emission, contains a smaller primary beam 
(5' as opposed to - 30' at 1384 MHz) and an RMS 
noise limit below the expected flux shown in Figure^ 
The image does not contain any significant continuum 
emission, and hence is not shown here. 

3.2 IRAS 15596-5301 

Figure [3] shows the total intensity emission at 2368 
MHz overlaid with 2368 MHz ATCA contours. There 
is a source coincident with the 4800 MHz source from 
[Caray et al.| (|2002 ) . The peak fluxes for this source 
are 19.6 and 14 mJy beam"^ at 1384 and 2368 MHz 
respectively Integrating within a region corresponding 
to the 5-cr level at each frequency, the total integrated 
fluxes for this source are 59, 53.7 and 49 mJy at 843, 
1384 and 2368 MHz respectively 

4 Discussion 
4.1 G333.125-0.562 

Figure |4] shows the expected synchrotron emission as- 
suming a local CR spectrum, for a variety of magnetic 
field strengths. It can clearly be seen that a field of 1 
mC would provide a 5-cr detection, whilst for a field 
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Figure 3: IRAS15596-5302: ATCA 2368 MHz to- 
tal intensity image overlaid with (black) 2385 MHz 
contours at 0.4 0.8 1.6 3.2 6.4 8.0 and 10.0 mJy 
beam~^ with the beam in the lower left-hand cor- 
ner. The same transfer function as Figure [2] image 
has been applied. 



strength of 500 /iC we predict a flux of ^-^0.5 mJy at 1 
GHz. Additionally, assuming no thermal contamina- 
tion, a 500 /xG field strength would produce a probable 
detection at 8640 MHz. 

Since there is no emission at 843, 1384, 2368 or 
8640 MHz down to the 1-a noise level of 0.5 mJy 
beam"^ at 1384 and 2368 MHz; 1 mJy beam"^ at 843 
MHz; and 50 /xJy beam~^, we place an upper limit on 
the magnetic field of 500 /iC. Despite being lower in 
frequency (hence having the expectation of a higher 
synchrotron emissivity), the upper limit on the mag- 
netic field is derived from the archival ATCA data, 
since the SUMSS catalogue has a limiting pe ak bright 
ness for sources 6 > —50° is 6 mJy beam~^ (Mauch et 
al.|2003t : an order of magnitude worse than the limits 



derived from our images and even worse still that the 
archival data at 8640 MHz. 



4.2 IRAS 15596-5301 

Figure [5] shows the full spectrum of IRAS15596. Two 
components are modelled: free-free emission, tracing 
the emission from the HII regions; and a modified 
Blackbody cur ve, following the dus t and infra-red emis- 
sion. Following [Garay et al. (2002), the free-free emis- 
sion was modelled using an emission measure (EM ^ 
Tied pc cm~^), for an electronic temperature of 10,000 
K. The higher frequency emission was fitted using two 
modified blackbody curves, of the form: 



5, = B.(Td)[l-exp(-T.)]Q. 



(12) 



where Bi^{Td) is the Planck function, Td is the dust 
temperature, and Qs is the solid angle subtended by 
the dust emitting region (see figure [s]) . The optical 
depth, Ti/, is allowed to vary with frequency; — 
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Figure 4: Predicted flux from the G333 cold core 
for different magnetic field strengths: 100 jaG 
(solid), 250 /iG (dot-dash), 500 fiG (dash), and 
1 mG (dotted). The crosses are (from left) the 
l-cr noise levels of 0.5 mJy, 0.5 mJy and ~ 50 
/iJy beam~^ derived from the 1384, 2368 and 8640 
MHz datum respectively. The physical parameters 
used in the determination of the flux densities are 
listed in Table 1. 



(ly/iyo)^, where is the frequency at which the op- 
tical depth is unity, and / 3 = 2, consistent with high 
mass star forming regions (Garay et al. 2002). The pa- 



rameters used to model the dust emission can be found 



Garay et al. (2002) 



From the integrated fluxes presented in the previ- 
ous section, it can be seen that the radio continuum 
emission traces the optically thin, free-free emission 
well with a spectral index between all radio contin- 
uum frequencies (843 to 4800 MHz) of a - -0.1. We 
have placed upper limits on the magnetic field in this 
cloud by turning the original question around, that is: 
how high can the magnetic field be without significantly 
altering the optically thin free-free emission spectrum 
that is observed?. Placing a conservative 10% charac- 
teristic error on the fluxes (due to RMS noise of the im- 
age, plus systematical uncertainties), a magnetic field 
is implied to be of the order of 0.5 mG is required to 
produce ^5-6 mJy of flux at 1 GHz. Taking 0.5 mG 
as the upper limit, t his is not incon sistent with the 
scaling relation from ( Crutcher|1999 ). 



5 Conclusion 

We undertook a search for synchrotron emission from 
molecular clouds due to production of secondary e^. 
We have modelled the expected flux using the local cos- 
mic ray spectrum, and shown that the emission should 
be observable above the more smoothly distributed, 
diffuse synchrotron emission of the galaxy if the mag- 
netic fields in the clouds are sufficiently high. We had 
initially selected three clouds, Sgr B2 (Jones, et a/., in 
preparation) because of the expected increase in cos- 
mic ray density near the GC, and two other clouds 
because they are isolated, dense, cold, dark clouds. 



Figure 5: Spectral Energy Distribution (SED) for 
IRAS 15596. The data points are: stars repre- 
sent the 843, 1384, 2368 and 4800 MHz SUMSS 
and ATCA data (integrated as described in the 
text); triangle is the 1.2 mm emission from ([Garay 
et al. |2002); the squares are the IRAS fluxes as re- 
ported in ( Garay et al.|2002 ). The fluxes modelled 
are the free- free emission (soli d curve) using emis - 
sion measures as described in ( Garay et al.||2002 ), 
and electron temperature of 10,000 K. The mod- 
ified blackbody curves are described in §5.2 and 
are for dust temperatures of 27 (dash) an d 100 K 
(dot-d ash), using the parameter set in ^ Gar ay et al. 



( i2QQ2| . 



with little evidence of star formation, and little or no 
previously detected radio emission. Additionally, their 
smaller distance makes the assumption of a local CR 
spectrum safer. For these nearby clouds, our conclu- 
sions can be summarised as follows: 



1. 



G333. 125-0.562, from 1.2mm, MSX and IRAS 
observations, was designated as a starless core 
( Garay et al.|2004 ). We have presented the first 
high resolution, radio continuum observations of 
this core, and do not find any radio emission 
at any frequency. This is in line with previous 
determinations that this cloud is not ob served 
at any wavelength longer than 1.2 mm (Lo et 
aLl|2007|. We place a 1-a upper limit on the 



2. 



flux from this cloud of 500 /xJy beam ^ at 1384 
and 2368 MHz, and 50 fiJy beam"^ at 8640 
MHz. The flux density upper limits in turn, were 
used to place an upper limit on the magnetic 
field of this cloud of 500 /iG, consistent with 
the magnetic fie ld - density correlation found by 
( Crutcher||1999 ). 

From 4800 MHz ob servations taken with ATCA, 
( Garay et al.]2QQ2| ) concluded IRAS 15596-5301 
to be a molecular cloud undergoing massive star 
formation. The flux at lower frequencies (includ- 
ing the 843 MHZ SUMSS data), presented here, 
confirms this. The spectral index is consistent 
with free-free emission from an HII region - in- 
dicating that this cloud is in an advanced stage 
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of star formation. Additionally, we constrain the 
magnetic field to be below 0.5 mG, a value which 
is also consistent with the scaling relation found 



in Crutcher ( 1999 ) 



3. Though we failed to detect synchrotron emis- 
sion from either cloud, we note that the initiat- 
ing particles are of ^-^GeV energies. If there are 
sufficient numbers, as our production spectrum 
indicates, then the upcoming GLAST mission 
may be able to probe the distribution of the CR 
spectum in our galaxy by selecting a sample of 
suitable clouds to search for 7-ray emission from 
7r° decay. 

We conclude that we did not find any evidence 
for synchrotron emission due to in these clouds. 
We note, however, that thermal emission from already 
formed massive stars, and nearby thermal emission 
tends to confuse possible detections, due to the low 
levels expected from the synchrotron emission. 
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